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S2CTI0N  1 


INTRODUCTION 

This  report  describes  the  problems  that  have  been  selected  to 
demonstrate  the  modeling  capabilities  of  the  SINGER  computer  code. 

The  problems  have  been  chosen  to  show  that  the  computer  code  has  the 
ability  to  predict  the  large  deflection,  finite  strain,  static  and 
dynamic  response  of  two  dimensional  frames  and  rings  composed  of 
straight  line  beam-column  elements. 

The  problems  have  been  selected  so  that  each  has  an  independent 
solution  check.  This  limits  the  scope  of  the  Individual  problems 
(especially  when  the  independent  check  Is  a  hand  calculation)  and 
each  problem  Is  able  to  check  only  a  few  features  of  the  program. 
However,  as  a  group,  the  problems  permit  checking  all  major  features 
of  the  program  and  have  proven  effective  In  developing  a  program 
which  can  model  relatively  complex  Interacting  inelastic  and  nonlinear 
effects. 

The  next  section  briefly  describes  each  of  the  four  problems 
selected  and  gives  a  comparison  of  the  results  obtained  by  SINGER 
with  those  obtained  by  the  Independent  check.  The  data  sets  used  in 
executing  these  problems  are  given  In  an  appendix  to  the  report. 


5 


SECTION  2 


DEMONSTRATION  PROBLEMS 


The  four  problems  selected  are:  (1)  plastic  beam,  (2)  elastic 
ring, (3)  elastics,  and  (4)  dynamic  rod.  Table  1  indicates  the  program 
features  tested  by  each  problem  and  the  source  of  the  independent 
solution  checks  that  will  be  used  to  evaluate  the  accuracy  of  the 
SINGER  analysis. 


2.1  PLASTIC  BEAM 

The  steel  wide  flange  beam  of  Figure  2.1  is  loaded  with  a  con¬ 
centrated  load,  P,  in  the  direction  of  the  y-axis  on  the  free  end. 

The  loading  is  monotonically  increased  so  that  system  deformations 
and  strains  pass  through  the  elastic  range  into  the  inelastic  range. 

The  beam  is  cantilevered  and  has  a  length,  L,  of  144  inches.  It 
has  a  flange  width,  b,  of  6  inches;  flange  thickness,  tf,  of  1  inch; 

web  thickness,  t  ,  of  1  inch;  and  a  total  depth,  h,  of  12  inches, 
w 

This  beam  was  analysed  with  two  different  stress-strain  curves  as  shown 
in  Figure  2.2.  In  each  case,  the  yield  point  is  33,000  psi,  and  Young' 
modulus  is  29,000,000  psi. 

In  the  elastic  range,  vertical  tip  deflections  were  calculated 

according  to  the  relation 

3 


tip 


PL" 

3EI 


(2.1) 


where  E  is  Young's  modulus,  and  I  is  the  moment  of  inertia  about  the 

4 

z-axis  which  for  this  beam  is  447.33  in. 

In  the  inelastic  range,  an  analytical  solution  for  deflections 
was  made  for  the  elastic-perf ectly  plastic  material  (Figure  2.2a)  and 
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TABLE  1 


DEMONSTRATION  PROBLEMS 


Number  Name 


Features  Tested 


Check 


1 

Plastic 

Nonlinear  material  model. 

exact 

Beam 

static  analysis  equation 
solution,  stress  analysis 

solution 

2 

Elastic 

Coordinate  transformations, 

exact 

Ring 

symsMtric  displacement 
constraints,  distributed 
loadings,  multldagree  of 
freedom  capability 

solution 

3 

Elastics 

Large  deflections,  adaptlva 
transformations,  geometric 
modeling  tests 

numerical 

solution 

(Huddleston) 

4 

Dynamic 

Cyclic  loading,  failure  model. 

exact 

Rod 

dynamic  analysis  capability, 
stress  modeling  limits. 

solution 
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for  the  material  in  which  strain  hardening  was  allowed  (Figure  2.2b). 

For  a  given  load,  the  moment  diagram  is  known.  The  stress  block  at  the 
support  is  determined  by  equating  the  integral  of  the  statical  moment 
of  the  stresses  to  the  applied  moment.  The  height  of  the  elastic  portion, 
y^,  is  then  determined  and,  from  similar  triangles,  the  strain  at  the 
extreme  fiber  is  found.  The  curvature,  'l' ,  is  obtained  by  dividing  the 
extreme  fiber  strain  by  the  depth  of  the  neutral  axis.  The  tip  deflection 
is,  by  moment-area  principles,  equal  to  the  moment  of  the  curvature 
diagram  about  the  free  end  of  the  cantilever.  It  should  be  noted  that 
this  analysis  does  not  consider  the  effect  of  axial  deformation. 

The  yield  load  was  found  to  be  17,086  lbs,  and  the  collapse  load 
for  the  elastic-perfectly  plastic  case  was  determined  as  20,850  lbs. 

This  collapse  load  was  obtained  by  setting  y^  of  Figure  2.3  to  zero. 

In  the  analysis  with  an  elastic-perfectly  plastic  stress-strain 
curve,  the  beam  was  divided  into  several  different  elements  of  varying 
lengths.  Studies  of  the  beam  as  1,  2,  3,  and  6  elements  with  three  load 
steps  were  made.  The  step  loadings  were 82%,  99%,  and  117%  of  the  yield 
load. 

The  results  were  more  accurate  in  the  cases  in  which  the  beams 
were  divided  so  that  the  energy  participation  of  the  elements  was 
evenly  distributed.  Output  from  SINGER  includes  this  information  on 
energy  participation.  Elements  with  a  high  percentage  of  participation 
are  simply  subdivided.  The  6  element  beam  gave  excellent  results,  but 
a  2  element  beam  with  a  114  inch  member  and  a  30  inch  member  at  the 
support,  where  yielding  began,  gave  comparable  results.  In  the  beams 
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that  the  elements  were  not  selected  to  balance  the  energy  distribution 
and  also  the  beam  considered  as  one  element,  it  was  found  that 
yielding  would  occur  during  the  99%  load  step. 

In  analysing  the  beam  with  strain  hardening  allowed  as  shown  in 
Figure  2.2b,  eight  elements  were  selected.  The  elements  were  arranged 
with  the  smaller  ones  at  the  support  and  the  longer  ones  at  the 
free  end  so  that  the  energy  participation  of  the  elements  would  be 
balanced.  The  SINGER  results  are  shown  in  Table  2  and  in  Figure  2.4 
as  load  versus  deflection  with  non-dimensional  units . 

The  SINGER  results  of  the  inelastic  beam  were  in  close  agreement 
with  the  analytical  results,  especially  if  the  elements  were  chosen 
so  that  the  energy  participation  was  balanced.  It  was  seen  that  a 
system  of  few  elements  properly  chosen  will  give  results  comparable 
to  a  system  with  a  large  number  of  elements.  When  strain  hardening  is 
allowed,  the  SINGER  results  vary  from  the  analytical  results  as  shown 
in  Figure  2.4.  The  SINGER  results  differ  from  the  analytical  results 
since  the  analytical  results  do  not  Include  the  effect  of  axial  de¬ 
formation. 

2.2  ELASTIC  RING 

The  steel  ring  of  Figure  2.5  is  symmetrically  loaded  by  a  uniformly 
distributed  internal  pressure,  p.  The  loading  is  chosen  so  that  system 
deformations  and  strains  that  are  within  the  elastic  range  and  linear 
analysis  assumptions  are  valid. 

The  ring  has  a  radius,  r,  of  36  inches;  thickness,  t,  of  2  inches; 
and  a  width,  b,  of  1  inch.  The  material  is  linearly  elastic  and  has  a 
Young's  modulus  of  29,000,000  psi.  The  internal  pressure  is  taken  as  2000 
psl  to  keep  the  stresses  well  within  the  elastic  range. 
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TABLE  2 


CANTILEVERED  BEAM  WITH  A  VERTICAL  LOAD  ON  THE  FREE  END 

P  -  17,086  lbs  L  -  144  in. 

y 


P/P 

y 

Tip  Deflection,  Inches 

Analytical* 

SINGER 

0.82 

1.076 

1.076 

1.00 

1.313 

1.313 

1.22 

5.310 

3.900 

1.67 

33.337 

32.162 

1.94 

61.516 

48.002 

2.07 

79.976 

54.636 

♦neglects  effect  of  axial  deformation 


Since  the  ratio  of  the  ring  thickness  to  its  radius  is  small,  uniform 


normal  stresses  can  be  assumed,  and  thereby  no  bending  moments  exist  in 
the  ring.  Also,  because  of  symmetry  in  geometry  and  loading,  no  shear 
forces  are  present.  From  the  well  known  ring  equations,  the  hoop  force 
is  given  by 


(2.2) 


and  the  radial  displacement  is 


u 


exact 


2 

Et 


(2.3) 


where  E  is  Young's  modulus.  For  internal  pressure  of  2000  psi,  the 
hoop  force  for  the  example  problem  is  72,000  pounds  and  the  radial  dis¬ 
placement  is  0.04469  inch. 

In  idealizing  the  ring  for  SINGER,  a  number  of  straight  line  segments 
in  the  first  quadrant  were  used  as  shown  in  Figure  2.6.  The  ring  quadrant 
was  restrained  from  rotation  and  axial  displacement  at  the  boundaries  to 
simulate  a  complete  ring.  A  unifomly  distributed  load  P  of  2000  pounds/inch 
was  imposed  normal  to  the  axis  of  each  segment  to  represent  internal 
pressure. 

To  determine  the  number  of  segments  required  to  adequately  represent 
the  ring  section,  data  sets  with  the  quadrant  divided  into  3,  4,  6,  9,  and 
18  segments  were  analysed  by  SINGER.  The  same  geometric  and  material 
properties  were  used  for  each  data  set.  The  only  variation  among  the  data 
sets  was  the  number  of  segments. 

The  results  of  the  analyses  are  shown  in  Table  3.  The  hoop  force,  T, 
and  the  radial  displacement,  u,  were  in  close  agreement  with  the  actual 
values  calculated.  Since  the  ring  sections  were  divided  into  beam  sections, 
SINGER  gave  some  shear  and  bending  moment  at  each  Joint,  but  these  values 
were  negligible  (five  orders  of  magnitude  smaller  than  the  axial  force) 
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The  accuracy  of  the  SINGER  results  as  compared  with  the  actual  calculations 


Is  shown  In  Figure  2.7.  In  the  figure,  F  and  F  can  be  either  T  and 

GXaCk 

T  ^  or  u  and  u 
exact  exact 

Since  the  results  for  all  the  data  sets  have  a  relative  error  of  less 
than  3.5%,  it  can  be  stated  that  SINGER  models  coordinate  transformations, 
displacement  constraints,  and  distributed  loadings  well.  The  analysis  of  the 
18  segment  ring  section  shows  that  SINGER  can  also  handle  multi-degree  of 
freedom  systems.  It  can  be  concluded  that  for  analyzing  rings,  It  Is 
sufficient  to  model  a  quadrant  with  9  straight  line  segments.  It  can  be 
seen  In  Figure  2.7  that  this  gives  results  comparable  to  those  with  twice 
as  many  segments  and  with  considerably  less  input  and  computing  effort. 


2.3  ELASTICA 


The  simply  supported  rod  of  Figure  2.8  is  subjected  to  the  slowly 
applied  axial  load,  F,  and  end  moments,  M.  The  load  is  Increased  so  that 
large  deflections  occur.  Since  the  load  is  near  the  Euler  buckling  load, 
the  beam-column  will  have  nonlinear  geometrical  behavior. 

The  rod  In  Figure  2.8  has  a  length,  L,  of  15  feet  and  cross  sectional 
dimensions  of  0.2  Inches  by  0.2  inches.  The  Young's  modulus  for  the  material 
Is  29,000,000  psi,  and  there  Is  no  yield  stress  so  that  the  material  is 
perfectly  elastic. 

The  "elastica"  problem  was  proposed  by  Euler  in  1744.  A  closed  form 
solution  for  this  problem  cannot  be  obtained,  but  Huddleston*  provides  an 
Independent  numerical  solution.  The  classical  approach  to  solving  the 
"felastica"  problem  is  to  solve  a  nonlinear  differential  equation  in  terms 

^Huddleston,  J.  V.,  "A  Numerical  Technique  for  Elastica  Problems."  Journ. 

Struct.  Div.  ASCE,  Vol.  EM5,  Oct.  1968,  pp.  1159-1165. 
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0.04 


Figure  2.8  -  Thin  Elastic  Rod  Subjected  to  Large  Axial  Force  and  End  Moments 


of  elliptic  integrals.  Huddleston  uses  a  numerical  method  in  solving  the 
nonlinear  boundary-value  problem.  In  particular,  he  employs  a  "shooting" 
method,  where  initial  conditions  are  assumed,  integration  is  performed  by 
a  predictor-corrector  method,  and  terminal  conditions  are  computed.  An 
iterative  process  is  used  to  adjust  the  initial  conditions  until  the 
boundary  conditions  are  satisfied . 

Figure  2.9  shows  the  results  of  Huddleston's  solution  for  a  thin  rod 
loaded  with  an  eccentricity,  e/L,  of  0.01.  Dimensionless  load  is  plotted 
against  dimensionless  deflection  of  the  center  line.  To  non-dimensionallze 
the  load,  the  Euler  buckling  load,  Pg,  was  used.  The  longitudinal 
deflection  at  the  roller  versus  the  axial  load  is  shown  in  Figure  2.10. 

The  SINGER  analysis  results  were  also  non-dimensionallzed.  The 
loading  was  stepped  in  11  increments  of  the  Euler  buckling  load  reaching 
a  maximum  of  94 1.  These  results  for  center  line  deflection  and 
longitudinal  deflection  are  shown  in  Figures  2.9  and  2.10  and  in  Table  4. 

In  performing  the  SINGER  analysis,  examples  with  the  rod  divided  into 
different  numbers  of  elements  were  studied.  It  was  found  that  the  best 
results  were  obtained  from  the  examples  which  contained  a  large  number  of 
elements,  in  particular,  those  with  30  elements.  The  results  were  better 
in  the  examples  where  the  elements  were  smaller  at  the  center  of  the  rod 
and  larger  on  the  ends  of  the  rod.  This  variation  in  element  lengths 
distributed  the  energy  participation  more  evenly  than  in  the  examples  with 
equal  element  lengths. 

The  comparison  of  the  Huddleston  and  SINGER  results  in  Table  4  shows 
that  SINGER  is  capable  of  handling  large  deflections.  It  can  be  concluded 
that  the  SINGER  code  can  simulate  systems  whose  equilibrium  position  depends 
on  adequate  representation  of  geometric  non-linearities. 


Figure  2.10  -  Non-Dimensional  Comparison  of  Load  and 
Longitudinal  Deflection. 


TABLE  4 


DEFLECTION  RESULTS  FOR  THE  ELASTICA  PROBLEM 


UNITS :  INCHES 


P/PE 

Center- 

SINGER 

-line  Deflection 

Huddleston 

Longitudinal  Deflection 

SINGER  Huddleston 

0.1 

0.25 

0.27 

0.01 

0.18 

0.2 

0.56 

0.63 

0.01 

0.54 

0.3 

0.96 

0.90 

0.01 

0.90 

0.4 

1.50 

1.26 

0.03 

1.26 

0.5 

2.25 

2.16 

0.07 

1.44 

0.6 

3.39 

3.33 

0.16 

1.62 

0.7 

5.27 

5.40 

0.39 

1.80 

0.8 

8.96 

9.00 

1.12 

3.60 

0.9 

18.17 

18.00 

4.64 

7.20 
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2.  A  DYNAMIC  ROD 


Figure  2.11  depicts  a  uniform  mass-less  rod  attached  to  a  support 
with  a  concentrated  mass  at  the  free  end.  The  constitutive  law  of  the 
material  is  defined  in  Figure  2.12.  At  time  zero,  the  unstrained  state 
of  the  rod  (x-o)  is  excited  by  an  impluse,  which  is  large  enough  to 
induce  inelastic  deformations  but  does  not  cause  fracture.  Consequently, 
the  rod  undergoes  cyclic  motion,  dissipates  energy  during  inelastic 
response,  and  attains  harmonic  motion  as  a  limit  state. 

The  purpose  of  this  problem  is  to  demonstrate  the  ability  of  SINGER 
to  predict  inelastic  cyclic  response  of  systems.  Since  the  response  of 
the  rod  to  an  initial  impulse  can  be  expressed  in  explicit  form,  it  serves 
as  a  basis  for  assessing  the  quality  of  the  response  predictions  produced 

by  SINGER. 

The  internal  force-displacement  relation  of  the  rod  is  piecewise 
linear.  The  corresponding  displacement  functions  can  be  described  by  two 
general  expressions.  If  the  stiffness  of  the  rod  is  nonzero  during  the  i 
linear  internal  force-displacement  domain  (Figure  2.13),  the  displacement 
function  assumes  the  form 

v^ 

x(t)  “  -d1  +  (xt  +  d±)  cosu^x  +  —  sinujjT  ,  (2. A) 


where 


d  -  7-  - 


i  k 


i 

AE 


i 

T  ’ 


wi  * 


T  -t  -  t^, 


(2.5) 


(2.6) 

(2.7) 


(2.8) 
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m  ■  10  lb  sec  /in. 
A  -  0.001  in.2 
£  ■  2880  in. 

<Jy  ■  36,00  psi 


Figure  2.11.  Dynamic  Rod 


Vol.  3,  User's  Guide 
for  further  details. 


Figure  2.12.  Stress-Strain  Curve 


Figure  2.13.  i  Response  Domain 
(x.  <  x  <  x,^) 


Figure  2.14.  Stress-Strain  Response 


and 


(2.9) 


In  Equations 2. 4-2. 8,  x  denotes  the  displacement  of  the  mass  m  from  the 

unstrained  state;  x^  and  x^+^  denote  the  displacements  at  the  beginning 

and  end  of  the  ifc^  linear  response  domain;  r^  is  the  internal  force 

acting  on  the  mass  at  the  beginning  of  the  1^  response  domain;  E^,  k^, 

represent  the  slope  of  the  stress-strain  curve,  the  stiffness  of  the 

rod,  and  the  natural  circular  frequency  during  the  I*1*1  response  domain, 

respectively;  k,t  denote  the  area,  length  of  the  tod,  respectively; 

t  is  the  time  measured  from  the  applcation  of  the  impluslve  excitation; 

t^  and  t^+^  denote  the  times  at  the  beginning  and  end  of  the  1th  response 

domain,  respectively;  and  x  denotes  time  measured  from  t^.  If  the  stiffness 

th 

of  the  rod  is  zero  during  the  1  response  domain,  the  displacement  function 
is  described  by 

x(x)  -  xA  +  v^x  -  y  ~  t2  (2.10) 

th 


where  v^  denotes  the  velocity  of  the  mass  m  at  the  beginning  of  the  i 
response  domain. 

On  the  basis  of  Equation  2.4  or  2.10,  other  measures  of  response 
can  be  obtained  as  follows: 

e(x)  -  ,  (2.11) 

o(x)  is  determined  by  the  strain  history,  Equation  2.11, and  the  stress- 
strain  relation  in  Figure  2.12,  and 

t(t)  ■  A  o(x)  .  (2.12) 

The  exact  response  to  a  specific  impulsive  excitation  (x^  -  0, 

■  10  in/sec)  is  depicted  in  Table  5  and  Figure  2.14. 
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A  comparision  of  the  response  predictions  produced  by  SINGER  with  the 
exact  response  is  presented  in  Figures  2.15  and  2.16.  Figure  2.15 
depicts  the  equilibrium  path  of  the  rod,  and  Figure  2.16  describes  the 
stress  history  of  the  rod.  The  agreement  between  the  exact  response 
and  the  response  prediction  of  SINGER  is  excellent. 
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Figure  2.16.  Stress  History  of  Rod 
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